Simple Cellular Model of 
Long- Range Multiplicity and pt Correlations 
in High-Energy Nuclear Collisions 

V.V. Vechernin and R.S. Kolevatov 

High Energy Physics Dep., St. Petersburg State University, 
1 Ulianovskaya str., 198504 St. Petersburg, Russia 



Abstract 

A simple cellular model for the description of the long-range multiplicity and 
Pt correlations in high-energy nuclear collisions originating from the string fusion 
model is proposed. Three versions of the model: without fusion, with local and 
with global string fusion are formulated. 

A Gauss approximation which enables explicit analytical calculations of the 
correlation functions in some asymptotic cases in the framework of the model is 
developed. The assumptions of the model and the validity of a Gauss approxi- 
mation are checked up in the simplest (no fusion) case when the explicit solution 
of the model can be found. 

The role of the size of cells is anylised. The modification of the results in the 
case of non-Poissonian distributions is also discussed. 







1 Introduction. 



The colour strings approach [jT], ^ is widely applied for the description of the soft part 
of the hadronic and nuclear interactions at high energies. 

In the frame work of this approach the string fusion model was suggested 0] . Later 
it was developed and applied for the description of the long-range multiplicity 

and pt correlations in relativistic nuclear collisions 0-0. 

The aim of the present paper is to formulate some simple cellular analog of the 
model, which enables explicit analytical calculations of the correlation functions in 
some asymptotic cases and drastically simplifies calculations in the case of real nucleus 
collisions. 

We check up the assumptions of the cellular model and the validity of a suggested 
Gauss approximation in the simplest (no fusion) case when the explicit solution of the 
model can be found. 

The paper organized as follows. Next section is devoted to the formulation of the a 
cellular analog of the string fusion model. The version of the model with a local string 
fusion is considered. 

The section 3 deals with the no string fusion limit of the model. The correspondence 
with the previous results JlOl is demonstrated. 

In the section 4 a Gauss approximation for the correlation function calculations is 
formulated. The results of the calculations in this approximation is compared with 
exact solution, which can be found in the no fusion case. 

In the section 5 the role of the size of cells is anylised. The cluster size dependence 
is considered and the version of the model with global string fusion is formulated. 

The modification of the results in the case of non-Poissonian distributions is dis- 
cussed in the section 6. 



2 Cellular approach to the string fusion phenome- 
non. 

Let us consider the collision of nuclei in two stage scenario when at first stage the 
colour strings are formed, and at the second stage these strings (or some other (higher 
colour) strings formed due to fusion of primary strings) are decaying, emitting observed 
particles. 

We'll consider three possibilities: without string fusion, with local and with global 
string fusion. The case with a local fusion corresponds to the model, where colour 
fields are summing up only locally and the global fusion case corresponds to the model, 
where colour fields are summing up globally - all over the cluster area - into one average 
colour field, the last case corresponds to the summing of the sources colour charges. 
(In section 5 we are refering to these cases as A) and B) correspondingly.) 

In the transverse plane depending on the impact parameter b we have some inter- 
action area S{b). Let us split this area on the cells of order of the transverse string 
size. Then we'll have M = S{b)/ao cells, where ctq = tttq is the transverse square of 
the string and tq ^ 0.2/m is the string radius. 

Local string fusion. At first let us consider the case with a local fusion. In this 
case the assumption of the model is that if the number of strings belonging to the z-th 
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cell is rji, then they form higher colour string, which emits in average jJ^oy/Vi particles 
with mean equal to p'^^/rji, compared with fiQ particles with (p^) = emitting by a 
single string. (Note that situation when some T^j = is also admitted.) 

Let us denote by nj and - the number and the average number of particles emitted 
by the higher string from i-th cell in a given rapidity interval, then 

ni = HoVVi (1) 

From event to event the number of strings rji in i-th cell will fluctuate around some 
average value - r/^. Clear that in the case of real nuclear collisions these average values 
Tjj will be different for different cells. They will depend on the position (s) of the i- 
th cell in the interaction area (s is two dimensional vector in transverse plane). To 
get the physical answer we have to sum the contributions from different cells, which 
corresponds to integration over s in transverse plane. 

The average local density of primary strings f}^ in the point s of transverse plane 
is uniquely determined by the distributions of nuclear densities and the value of the 
impact parameter - b. They can be calculated, for example, in Glauber approximation. 
We'll do this later in a separate paper. In present paper we consider that all rj^ are 
already fixed from these considerations at given value of the impact parameter - b. 

If we introduce: 

M M 

N = j2m, N = j2m (2) 

1=1 i=i 

then clear that N is the number of strings in the given event and is the mean number 
of strings for this type of events (at the fixed impact parameter b. 

To go to long-range rapidity correlations let us consider two rapidity windows F (for- 
ward) and B (backward). Each event corresponds to a certain configuration {rji, 
of strings and certain numbers of charged particles {rii, n^^^} emitted by these strings 
in the forward rapidity window. Then the total number of particles produced in the 
forward rapidity window will be equal to np: 

M 

nF = ^ni (3) 

i=l 

The probability to detect np particles in the forward rapidity window for a given 
configuration {q^, ...,rj^} of strings is equal to 

M 

where p^_{ni) is the probability of the emission of rii particles by the string rji in the 
forward rapidity window. By our assumption (|ip 

oo 

= E ^iPm^^i) = /^ov^ (5) 

ni=0 

If we denote else by W{rj^, ■■■^Vm) probability of realization of the string con- 
figuration {rj^, ...,rj^j} in the given event, then the average value of some quantity O 
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under condition of the production of np particles in the forward window will be equal 
to 

One has to omit in this M-fold sums one term, when all rji = 0, which corresponds to 
the absence of inelastic interaction between the nucleons of the colliding nuclei. 

If the O in the number of particles produced in the backward rapidity window ub 
in the given event, then (for (n^)^^ correlations) we have to use 

M 

i=l 

If the O in the mean squared transverse momentum of particles produced in the back- 
ward rapidity window p"^^ in the given event, then (for (Pt^)^ correlations) we have 
to use: 



i=i l^i=i y'li l^i=i y 'H 

Later we'll assume that numbers of primary strings in each cell rji fluctuate inde- 
pendently around some average quantities rj^ uniquely determined by the distributions 
of nuclear densities and the value of the impact parameter - h (see above), then 

M M 
i=l i=l 

For clearness we'll sometimes address to a simple "homogeneous" case, when all 
?7j (but not the r^j, which fluctuate!) is equal each other in the interaction area rj^ = 
Tj. The parameter rj coincides in this case with the parameter rj used in the papers 
IP, 1^, 1^ and has the meaning of the mean number of strings per area of one string 
(?7 = (mean string density) x do). In general case the parameters fj^ have the same 
meaning, but with mean string density depending on the point s in the transverse 
interaction plane {fj^ = (mean string density in the point s) x ctq)- 

If we assume else the Poissonian form of p^^{ni) {Pa{x) is the Poisson distribution 
with X = a): 

i'J-^) = P.oV^(-^.) - e--^^^^^ (10) 
then we find for -P{^^^...^^^}(^ir): 

P{,,,...,v,,}inF) = P,oE.^^np) (11) 

3 No string fusion. Correspondence with the pre- 
vious results. 

No string fusion. In the no fusion case we have the same formulae, but instead of 
(|), (|^ and (|^), we have to use 

oo 
ni=0 
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M 



M{r,„...,n,,},n^= f^oY.Vi (13) 



i=l 



4=1 2^i=l 'li 

which immediately leads to the absence of (Ptfi)™^ correlations. 
In this case we have to use also 

PvS^^) = P.ovM.) ^ e-'^^^^^J^ (15) 

and 

= P,oE.ri.M (16) 

instead of (|To|) and (11). Then for {n^)^^ correlations we find 
Introducing under the sums 

l = E'^^,E... (18) 

and putting everywhere = ^ '^^ 

, . _ f^oENNW{N)p^^^{n,) 

where 

WiN)= '^^,E,..n^(^0 (20) 

If we also admit the Poissonian form for w{rii) and that all fj^ = rf (the homogeneous 
case): 

= P,{m) = 6-"^^ (21) 

then we find 

W{N)=p^j^{N) (22) 

and ([T9|) has the form 

, \ _ /^O EjV A^PM^(A^)P;.oiv(^F) 
EivPMr,(^)PM0iv(^F) 

As = X^Af A^W^(iV) = Mr] we see, that (|23|) coincides with the formula from the paper 



T0[] , where the notation n = po for the mean multiplicity from one string (emitter). 



The and A^ is the number and the mean number of strings (emitters) was used. 

, , _nENNp^{N)p^NM . . 

ENpj^{N)PjiN{nF) 
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4 Gauss approximation. 

Let us now evaluate (|^) at iV 3> 1 and fioN 3> 1. At these assumptions we can replace 
J2n by / dN and the Poissonian distributions by Gaussian distributions ((yf^^(a;) is the 
Gauss distribution with x = a and x"^ — x"^ = a"^) . 

(n„-(n„>„)2 



P^.oNM ^ 9{n^)^,aA''F) = ^^^ ^ (25) 



with o"!' = {np)j^ = hqN. Similarly 

-1 (Af-iV)^ 



M^) ^iv,. J^) = (26) 



with 



and 



■3- 

27ro" TV 

with (T^ = A^. 
So we find 

^'^B)n, = J^^_^g-^(7V,n,) (27) 

^(iV,M^(^^+(M^ (28) 

^ = £ - 1 + ^ - f29l 
dN N 2 2/ioAr2 ^ ^ 

Let us denote by A"* the point where ^ = 0. Then we can evaluate (^) as follows 

(^B)n^ = /^oA'*(«f) (30) 
In relative variables we can rewrite condition = as follows 

dN 

z'-z^ = ^{f-z') (31) 

which defines z as function of /, where z = N^/N and / = Up/ (up) = np/{fiQN) and 
then 

= fioN^irip) = fioNzif) = {np)z{f) (32) 
Let us define correlation coefficient as 

b=^^^\ -^1 (33) 



Prom (31) we have 

dz fiof 



df 3^2 + z{no- 2) 

and then 



Because as clear from (|3l|) at / = 1 one has z = 1. 



(34) 
(35) 
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In reality for one string /iq = ^^D-, where ^ ^ 1.0^ 1.2 and Ay is the width of the 
rapidity window. If one chooses backward and forward windows of a different width 
Ays 7^ A^/i?, then //qb 7^ /ioF, and instead of (^) we have 

rfn^ - /xoF + 1 ^ ^ 

or for "relative" quantities 



(37) 



We see that in the last case correlation coefficient depends only on value of /iof in the 
forward window (see physical explanation of this fact in the end of the section) . 

In Figs. 1-2 we present the results of the exact (dashed lines) and approximate (in the 
Gauss approximation) (solid lines) calculations of the function (n^)^^ using formulas 
(P^) and (^2p correspondingly at the different values of the mean number of strings 
A'" = Mr] and no in the no fusion case. We see that the Gauss approximation works 
very well starting from = 4 especially in the region rip = (rip), where the most of 
experimental points lay and where correlation coefficient is defined. 

One has also to pay attention on the approximate linearity of the correlation func- 
tions obtained here in the case of two Poisson distributions by use of the formula (p^, 
which coincides with the corresponding formula from UTOf, the problem dealt with in 
that paper. 

In Figs. 3 we present the results of the exact (dashed lines) and approximate (in 
the Gauss approximation) (solid lines) calculations of the correlation coefficient b = 

— -li^ — - |„ ^ as function of /iq using formulas (0) and (^) correspondingly at the 

different values of the mean number of strings = Mrf in the no fusion case. We see 
again that Gauss approximation works very well starting from = 4 for any /xq- 

As one can see from Figs. 1-3 we have A^-independence for correlation functions and 
correlation coefficient b starting very early (from A^ = 4). It's interesting to mention 
that as one can see from Figs. 4, 5 at that the resulting distributions P{np) are changing 
drastically with iV from iV = 4 to iV = 128. 

We see also that we have practically ideal Gauss distribution for P{np) at N = 128. 
This is in agreement with the experimental data. Unlike the case of pp-interactions 
(where NBD for P{np) takes place) in the case of nuclear collisions with a large number 
of emitting centers the ideal Gauss distribution for P{np) has been observed exper- 
imentally for central PbPb collisions (i.e. at fixed value of impact parameter b and 
hence fixed S{b)) (see, for example, Fig. 6 in |]lll). Note that for nuclear collisions in 



Glauber approximation we have also Gauss distribution for P{np) at fixed value of 
impact parameter b (see, for example, |T2| ). 

Physical interpretation. Let us to discuss in the conclusion of the section why the 
correlation coefficient for "relative" quantities b (|37| ) depends only on the multiplicity 
in the forward rapidity window. 

The correlations between and rip in the model under consideration arise only 
through fluctuations in the number of strings A^. At large values of /xqf ^ 1 it's more 
probable that fluctuation in the number of forward produced particles rip was caused 
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by the fluctuation in the number of strings N than by the fluctuations in rii and vice 
versa at small values of fioF <^ 1 it's more probable that fluctuation in the number 
of forward produced particles Up was caused not by the fluctuation in the number of 
strings but by the fluctuations in rij. 

Formally we can see it from analysis of the maximum A^.,, of the function ip{N,np) 
(^), in which the first term originates from the fluctuations of and the second term 
originates from the fluctuations of nj. If we have some fluctuation (i.e. rip ^ fJ'oN), 
then at /ioF S> 1 A^* — >■ rip/ fiQ and at /iof ^ N . 

So if every string emits in average large number of particles (/ioF ^ 1) in the given 
forward interval Ayp then based on information of rip we can do the justified conclusion 
on the number of strings A^ in the given event and hence expect corresponding change 
in the n^. And if every string emits in average small number of particles (/ioF ^ 1) in 
the given interval ^yp then based on information of rip we can't do any conclusions 
on the number of string A^ in the given event. 

Note that more detail analysis shows that this conclusion is not based on the specific 
(Poissonian) form of the distributions. At small n j <^ 1 one always will have CTj ^ Hi 
due to discrete nature of rii. 

5 Cell size, cluster size, global string fusion. 

Let us go back to the string fusion and consider the case, when Tc - correlation radius 
(cluster size) is not equal to the string radius - rg (7 = Tc/tq > 1), then the square of 
the cluster in a transverse plane will be equal to AS" = 7^(To and 7^ is the square of 
the cluster in string square units. If j-th cluster was formed in the given event by rrij 
strings, then it will have rjcj = rrijao/AS = rrij/'-f^, where j = 1, ...,Mc and Mc is the 
number of the clusters: Mc = S{b)/AS = 7^) = M/'j'^. (We keep notation M for the 
quantity S{h)/ao = M.) The S{h) is the transverse square of interaction area of nuclei 
at given impact parameter h. Note that in given approach we consider the clusters of 
the fixed (by hands) area (AS) with fluctuating from event to event number of strings 
rrij forming it (i.e. the rjcj fluctuates around r] in the homogeneous case). 

Now the elementary emitters will be not strings, but clusters. The mean number 
of particles emitted by such cluster will be equal to ^icy/fhj (in the case with string 
fusion), where Hc = jJ^o^ = l^ol"^ = ^^1/7^- So to study the dependence on the cluster 
size AS in the given model we have to 

1) increase the luminosity of elementary emitters /ioy^ ~^ f^cy/Vcj = /^o7^-\/'^j/7^ = 
/xo7^/m~ 

2) simultaneously reduce the number of clusters M — ^ Mc = S{b)/AS = S'(6)/((To7^) = 
M/72. 

Note that there is no sense to introduce the clusters AS > ctq (i.e. correlated 
fluctuations of the r]i within area AS*) in the case without string fusion. Similarly it's 
no physical reasons to consider clusters at small values of rj {rj < 1) even in the case 
with string fusion. 

On the contrary at large rj [r] >> 1) there are physical reasons in the case with 
string fusion to consider two possibilities: 

A) AS = (To, 7^ = AS/ao = 1, = M = S{h)/ao 

B) AS = S{b), 72 = AS/ao = M, = 1 (one cluster). 
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The first one corresponds to the model, where colour string fields are summing up only 
locally (formulated in section 2) and the second one corresponds to the model, where 
colour string fields are summing up all over the interaction area S{b) into one colour 
field. 

In the case A) at the first stage we have M = S{b)/crQ with rji, i = 1, ...,M fiuc- 
tuated around 77. Then we have to generate particles from each area o"o with average 
multiplicities equal to fio^/rji (see section 2). 

Global string fusion. In the case B) we must MODIFY our formulae, as at first 
stage we also have M = S{b)/ao (like in the case A)) with rji, i = I, ...,M fiuctuated 
around rj. Then (unlike the case A)) we have to find average y]c = j^HiVi = J4 given 
event, and then to generate particles from one cluster with average multiplicity equal 



to fic^/Vc = f^oMy/rj^ = fioM^N/M = ^qVMN. (Note that for configuration when all 
rji = r] the average multiplicity in the case A) and B) is the same yUo(>S'(6)/cro)y^. In 



general case they will be slightly different, as y jjJ2i J2i 
So in case B) we must replace @ by 

_ E(,,....,,,|(0)H,...,,^^|,.^H^(r/,, ...,v,,)p^^^^^^{n,) ^^^^ 

In the case B) we must ALSO MODIFY our expressions for ^m^^f ' ^^^^s 

of the backward production from configuration {rj^, ...,rij^j}. We have to use instead of 
(0) for correlations: 



with M = S{b)/ao and instead of (^ for {pIb)^,^ correlations: 



sS'ft (39) 



jfE'n (40) 



Again we see that the difference with the case A) consists in replace J2i \f^i ~^ 
\J Vi- As a consequence calculations in the case B) are much more simple as we 
can reduce all sums J2{-ri-^,...,rjj^j} to one sum J2n using identity (|I^) as in the no fusion 
case. 

So in the fusion case with one cluster (case B)) we can write simple formulas as in 
the no fusion case. Namely we have 

, , _ A.x/ME«v^^t^(jV)p„^^(»F) 

E„M'(iV)p„^^(,V) <"> 

and 
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with M = S{h)/aQ and W{N) is given by the formula (pO]). We see that in this case 
B) (one cluster at large r/) n-n and pt-n correlations are connected 

Note also that in the case A) for Poissonian distributions Pj^-{ni) ([l0| ) we have 
for -P{^^^...^^^}(^i7) the formula (pJ]) which is very similar to (p^) but with 



instead of y^J^iVi- More over we can fulfil in the case A) the summation (^) over 
{n^, n^} and get the similar formula for P^^^ rij^^yi^p) ^^^o in the cases of binomial 
and negative binomial distributions (see the next section). In these formulas the rate of 
Up production depends only on some average quantities for configuration {r^^^, ...,1]^^} 



(on J2i \fWi in the case A) and on \J -^Y^i Vi the case B)) So we expect that we'll 
have very similar results for n-n and pt-n correlations in the cases A) and B). 

6 Non-Poissonian ( binomial, negative binomial ) 
distributions. 



Important note: Even if we have no Poisson distributions ([T5|) for Prj^irii) we still 
can get Gauss type formula ( |25D for P|^^ nj^,j}i^F) at M 1 using (H) and the central 
limit theorem of the probability theory 

^ 9{n^),aA'^F) (44) 

with (up) = J^i^i and ap = J2i <^h ^ut now ap can be not equal to {up). 

Let us consider, as an example, the binomial distributions instead of Poisson dis- 
tributions. 

p,M^) = Pk.M) ^ ci^y^^ii - (45) 

Unlike Poisson distribution there are two parameters in binomial distribution. We 
choose the same parameter A for all p^. (rij) distributions. In this case we can sum 
these distributions like Poisson distributions due to formula 

M 

Pj:^k.,xM= E K,E.r^.I[Pk.^in^) (46) 

For binomial distributions we have 

Hi = kiX a1 = kiX{l - A) = ni{l - A) (47) 

In the no fusion case we also have = ^oVi and hence ki =ni/X =^rii. So we have 
to choose A so that to ensure the ^ will be integer, in other respects the parameter 
A is arbitrary. (Clear that n^"-^ = ki and the case of small A, when ki = n™"^ ^ n,, 
corresponds to the Poisson case.) Then we find 



= /5<vl Mf) (4^ 
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with {up) =Y,i'^i =/^o Y,i Vi '^F = J2i ^1 =(1 ~ '^)/^o Y,i Vi- We see that for binomial 
distribution a\ = {\ — X){np) with < A < 1. Note that case A ^ 1 corresponds to 
the situation, when all rij = n"^"-^ = ki and there are no dispersion all af = 0. 
Similarly for w{7]i) starting instead of (^I]) from 



HVi)=Pr^,xM) (49) 

with 

Vi = n\ = v = - K) = - \) (50) 

where now ?7/A^ must be integer. Then we find instead of (p^) 

W{N) = (3m, , (N) (51) 

with N =J2iVi =Mr] and =J2i<^^i =(1 — \)Mri and again = (1 — A^)A^ with 
< A^ < 1. Clear that the case A^ — » corresponds to the Poissonian case and the 
case A^ ^ 1 corresponds to the situation, when N = N and there are no dispersion 
a% = (fixed number of strings). 
Then instead of ( P^ ) we find 

^ iN)P^^^inp) ^^^^ 

T— ,A,) A ' 

V 

In Gauss approximation we have the same formulas (|25| ) and (^) but now with = 
{np)jy{l — A) = /ioA^(l — A) and = A^(l — A^). We have also the same formula (|27|) 
but now with 

Clear that the case A — > 1 corresponds to 6{np — ^qN) and = Up/fiQ. Each string 
emits exactly /io particles and the number of strings in the given event can be 
uniquely reconstructed on the value of rip which leads to 100% correlations (see dis- 
cussion in the end of the section 4). The case A,, — > 1 corresponds to 6{N — N) and 
N = N, which corresponds to the fixed number of strings. There are no dependence of 
on the value of Up in this case which leads to no correlations (note that we consider 
the no fusion case). 

We'll have all the same formulae (^,^,^) for (n^)^^ and b, but with modified 
equation for z{f): 

z'-z' = ^if-z') (54) 

where 

.^1^ (55) 

We see again that at A — 1 we have k — » oo and 2 = /; at A^ ^ 1 we have k — > and 
z = 1 (doesn't depend on /). Instead of (|34D and (|35| ) we have 



dz nonf 



df 3z^ + z{noK — 2) 



(56) 
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and then 

b = (57) 

We see that really at A 1 we have k, ^ oo and the correlation coefficient 6 = 1; at 
— i> 1 we have k — ^ and the correlation coefficient 6 = 0. 

Clear that we can do all the same for negative binomial distributions. We can also 
use any combinations of these distributions - one type for p^. (rij) and another type for 
w{rii). 



7 Conclusion. Acknowledgments. 

We see that suggested simple cellular model for the description of the pt and multiplic- 
ity correlations in high-energy nuclear collisions and the Gauss approximation which 
enables explicit analytical calculations of the correlation functions in some asymptotic 
cases give the adequate results in the no fusion case. 

In the next paper we plan to present the results of the correlation functions calcula- 
tions with taking into account the string fusion phenomenon, based both on numerical 
summations on the configurations {r]^, ...,rij^.j} using formulae of the section 2 and on 
analytical calculations using the Gauss approximation discussed in section 4. 

In conclusion the authors would like to thank M.A. Braun and G.A. Feofilov for 
numerous valuable discussions. The work has been partially supported by the Russian 
Foundation for Fundamental Research under Grant No. 01-02-17137-a. 
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muO=1.0, both Poisson, no fusion 
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Figure 1: {n^)^^ correlation functions in the no fusion case, A^-independence. The 
check of vahdity of the Gauss approximation. The {N) — N is the mean number of 
strings, the //q = A*of is the mean number of particles emitted by one string in the 
forward rapidity window, /iq — 1. (The p{ni) and w{'r)i) are both Poissonian.) 




Figure 2: The same as in Fig.l, but for /xq = 2. 
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correl.coeff., both Poisson, no fusion 
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Figure 3: {n^)^ correlation coefficient h (|37D in the no fusion function of /xq, 

A^-independence. The check of vahdity of the Gauss approximation. The (A^) = is 
the mean number of strings, the /iq = /ioF is the mean number of particles emitted by 
one string in the forward rapidity window. (The p{ni) and w{rii) are both Poissonian.) 




Figure 4: The same as in Fig.l, but plots for P{np) added (not normalized, arbitrary 
units) . 
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Figure 5: The same as in Fig. 2, but plots for P{np) added (not normalized, arbitrary 
units) . 
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